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Abstract. We study the formation of deeply bound pionic atoms in the (d,^He) reactions theoretically 
and show the energy spectra of the emitted ^He at finite angles, which are expected to be observed 
experimentally. We find that the different combinations of the pion-bound and neutron-hole states dominate 
the spectra at different scattering angles because of the matching condition of the reaction. We conclude 
that the observation of the (ci,'^He) reaction at finite angles will provide the systematic information of 
the pionic bound states in each nucleus and will help to develop the study of the pion properties and the 
partial restoration of chiral symmetry in nuclei. 
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1 Introduction 

Deeply bound pionic atoms were discovered in the (d,'^He) 
spectra in 1996 [Tl[51|3] following the theoretical predictions 
[HIS]. They have been considered to be one of the good 
systems to deduce the pion properties at finite density and 
to obtain precise information on the partial restoration of 
chiral symmetry in nuclei [S1[71IH] ■ 

So far, the {d,^He) spectra have been obtained in near 
recoilless kinematics to observe the peaks of the pionic 
state formation with a neutron hole in the quasi - sub- 
stitutional configurations [SJiniE]- The Pb and Sn isotopes 
were used as the target nuclei in the experiments [3l[6j|9] 
since the spectra were expected to show the peak struc- 
tures due to one dominant [tt ^ n~^] configuration based 
on the theoretical evaluations pilllll) . The simple struc- 
ture of the observed peak is important and effective to 
deduce the pion binding energies precisely from the ob- 
served spectra. Actually, in Ref. the binding energy 
and width of pionic atoms in Sn isotopes were determined 
precisely and the partial restoration of chiral symmetry 
was concluded based on the data. 

To develop the studies of pion properties and symme- 
try restoration in nuclei further, we think that we need to 
obtain improved information from experiments. Actually, 
the systematic errors coming from the uncertainties of the 
neutron distribution and the absolute energy calibration 
were reported in Ref. [5]. In addition, we need more sys- 
tematic information on bound states for unique determi- 
nation of the pion-nucleus interaction, which is required 
to fix the potential strength related to chiral symmetry. 
Recently, it was reported in Ref. [12] that the simulta- 
neous observation of various pionic bound states such as 
Is and 2s in the same nucleus may be helpful to reduce 



these errors and to develop our studies. Along with this 
line, in actual experiments, the pionic 2s state observation 
with the Is state has been expected in new high precision 
experiment in RIBF/RIKEN 12,13,14 . 

In this paper, we consider theoretically the new possi- 
bility to observe various pionic bound states in the same 
nuclei by observing the (d,'^He) spectra at finite angles 
together with the forward direction. These observations 
at finite angles are expected to be possible in the experi- 
ments at RIBF/RIKEN [T5]. As well-known in the nuclear 
inelastic scatterings, we can expect to have the manifesta- 
tion of different subcomponents of pion and neutron hole 
states at finite angles due to the matching condition with 
different momentum transfer, and expect to determine the 
binding energies and widths of various pionic states simul- 
taneously in each nucleus. The momentum transfer depen- 
dence of the effective number was studied in Ref. [5] for 
different incident energies. The angular dependence of the 
resonance peak structure will be helpful to determine the 
angular momentum quantum numbers of the populated 
states directly from the (c?,'^IIe) data. 



2 Formulation 

We modify slightly the theoretical model used in Refs. [5] 
[TUlfTT] to study the angular dependence of the (d,'^He) 
spectra. The modifications are implemented by including 
the kinematical correction factors K in Eq. (fTTj) as ex- 
plained below. A similar consideration was also given in 
Ref. [TB| for the {K~,N) reaction. 

We start our discussion by considering the d + n ^ 
■^Hc + TT^ reaction, which is the elementary pion produc- 
tion process for the formation of the pionic atoms in the 
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((i,^He) reaction. Our model is based on the relativistic 
cross section formula [17 with the Lorentz invariant pion 
production amplitude T (pd, Pm PHe, P-k) ■ The cross section 
can be written as, 

1 MdMn, , ,,2 
Urol rjn 



X (27r) 6 {pcl +Pn- PHc - Ptt) 



Mhc rfpHc 1 dp^ 



(1) 



where the subscript of the momentum p, the energy E 
and the mass M indicates each particle participating in 
the two body reaction d + n ^He + 7r~. The relative 
velocity Vj-d between d and n can be written as, 



1 



where s is the Mandelstam variable and A is the Kallen 
function. 

The elementary cross section in the center of mass 
(CM) frame and the laboratory frame at scattering angle 
^rfHc which is required for the effective number approach, 
are written as 



da 

dQ-nc 



CM 



I T (prf , p„ , pHc , Ptt ) P MdMn Mhq 



and 

/ da 



\dQ 



He 



lab 



clc 



\T{pd,P„ , PHc, Pvr) pMriA^n^Hc 



(3) 



EttIPhcI + -E^HcdPHcl - \Pd\ C0s6'<jHc) 



lab 



(4) 



where the superscripts 'lab' and 'CM' indicate the frame 
where the kinematical variables should be evaluated. The 
ratio of these cross sections can be written as. 



/ da \^'^^ , / da \ 

[dn^eJdJ [dn^oJd, 



CM 



E^IPHcI + -E^HedPHol ~ \Pd\ COsOdHc) 

AV2(,,Af2^,M^) 



lab 



This expression can be simplified as, 

lab / , \ CM 



/ da \^ , / da \ 



clc 



|„CM|2 • 
\PHo I 



(5) 



(6) 



at forward angles ddHc = 0° as shown in Ref. [4i[5j. The 
experimental elementary cross sections in the CM frame 
were reported in Ref. |18) . 



Then, we consider the pionic bound state formation 
cross section for the heavy nuclear target case in the same 
framework. In the effective number approach used in Refs. 
[51 ll0l[TT| . the elementary cross section is factorized from 
the bound state formation cross section as shown below. 
The bound state formation cross section for the heavy 
nuclear target case can be written in the laboratory frame 
as, 



V 



(27r) 



:dp 



He, 



(7) 



with the S matrix. 



Sfi — I dtdx 



Vv 

X iT{pd,Pn,p-ac,P7v) 




Xd{x) 



Mx), (8) 



where 0^ and ipn indicate the wavefunctions of the pion 
bound state in the daughter nucleus and the neutron bound 
state in the target nucleus, respectively. The wavefunc- 
tions of the projectile (d) and the ejectile (^He) are de- 
noted by Xhc ^-iid Xd- T{pd,Pn,PKc,P-!T) indicates the same 
meson production amplitude appeared in Eq. ([T|) and de- 
scribe the d + n '^He-|-7r~ transition. V and T indicate 
the spatial volume and the time interval for the transi- 
tion [17]. We assume that only the neutron in the target 
nucleus participates in the reaction and other nucleons are 
spectators. The formation cross section of the pion bound 
states in the laboratory frame can be expressed as, 



d^a 

di?Herf^He 



lab 



E 



1 MdM„Mn^ IphcI 
2(2^)2 E^E^ \pd\ 

r 1 



\T{pd,Pn,Plic,PTr)\' 



2n AE^ + n/4: 



dxxHcix)(t)lix)-)p„ix)xd{x) 



lab 



(9) 



where the subscript 'A' indicates the {d,^lie) reaction for 
the nuclear target. Here, we assume that the target nucleus 
is so heavy that we can ignore the recoil energy of the nu- 
cleus and that we can consider the CM and the laboratory 
frames to be the same. The (5-function responsible for the 
energy conservation is replaced by the Lorentz distribu- 

r 1 

tion function -—r: , to account for the width F 

27r AE^ + r2/4 

of the pion bound state, where AE is defined as AE — 
Ehc + Et^ — Ed — En. En is the energy of neutron in the tar- 
get nucleus and defined as En = Mn — Sn {jn ) with the neu- 
tron separation energy Sn from the j„ single particle level 
and the neutron mass M„. Et^ is the energy of the pion 
bound state defined as E.,^ = — with the pion 

mass and the binding energy B.E. of the bound state 
indicated by fjr- The reaction Q-value can be expressed as 
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Q = AE- + B.E.ii^) - SnUn) + (Af„ + - A/hc), 
where M„ + Md - Mho = 6.787 MeV. 

We express the effective number for each eigen state of 
the total spin of neutron-hole and pion-particle states in 
the daughter nuclei. We decompose the sum of the final 
states as 



Table 1. Elementary cross sections in CM frame {da /dS7)'^^ 
and those in lab frame (da/dfi)]^^ in unit of [/ib/sr]. The 
scattering angles 0][^^ and 6^^^ indicate the angles between 
deuteron and ^He in the d + n "^He + tt" reactions in CM 
and lab frames, respectively. We use Eq. (|14|l to fit the exper- 
imental data in Ref. [18] . 



E 



ph JM 



(10) 



where indicates the sum of all combinations of the pion 

ph 

bound states and neutron hole states in the final states, 



and 



JM 



indicates the sum of the states with different total 



angular momentum for each combination of bound pion 
and neutron hole state. By neglecting the residual inter- 
action effects for the final state energies, we can rewrite 
the bound state formation cross section as, 



lab 



/ da 



\rfi7Hc 



lab 



ph 



Ncs, (11) 



where we have used the elementary cross section 



/ da 



\ df2}ic 



instead of |Tp. The effective number A'off defined as. 



lab 



clc 



E 

JM 



drXHe{r)[(f>*iA'^) ® 'PjAr)]jMXdir) 



(12) 

We introduce the distortion effects to the wavefunctions 
Xjjg and Xd by Eikonal approximation as described in 
Refs. [Sl fTUllllj . The kinematical correction factor K is de- 
rived by factorizing the elementary cross section in the 
laboratory frame Eq. ^ from the bound state formation 
cross section for the heavy nuclear target case Eq. and 
is defined as. 



K 



IPhcI E,iEA 



Ehc IphoI - |Pd|cos6', 



lab 



(13) 



where A indicates the momentum and energy which should 
be evaluated in the kinematics of the nuclear target case. 
Since we consider the heavy target limit, the scattering 
angle of nuclear target case does not appear in Eq. ([T^ . 
It should be noted that in addition to the wavefunctions 
(j)^^ and i/ij,, , AE, F and K also depend on the final 'p/i' 
combination through the binding energy and width of the 
pion bound states and neutron separation energy. The K 
factor is evaluated by the energy and momentum at the 
kinematics satisfying the condition AE = for each 'p/i' 
combination. This correction factor is A' = 1 for the recoil- 
less kinematics at 9]^^^ = 0° with S'„ = and B.E. = 0. 

Here, we have assumed the same meson production 
amplitude T(pd,Pn,PHo,P7r) for the pion bound state for- 
mation process in the heavy nucleus as the elementary 



glab 



/da_y'"' 
\ dQ J clc 



0.0° 



1.0° 



2.0° 



3.0° 



,CM 
dHc 
da -CM 



\dnj,i, 



0.0° 
8.00 



15.0° 
7.48 



31.0° 
6.05 



49.5° 
4.02 



lab 



process by the impulse approximation. We need to know 
the microscopic mechanisms of the meson production to 
evaluate the possible modifications of the T amplitude be- 
yond the kinematical corrections. 

As a summary of this section, we use Eq. (TTT]) to calcu- 
late the angular dependence of the pionic atom formation 
spectra. This formula newly include the kinematical cor- 
rection factor K which has a certain contribution to the 
angular dependence of the cross section as we will see later. 

As for the elementary cross section, we obtain f — -- — 

Vc^^He/ele 

in Eq. ([TT|) by the formula Eq. ([5]) from the cross section 
in the CM frame which have been determined experimen- 
tally [H]. 



3 Numerical Results and Discussions 

First, we consider the angular dependence of the elemen- 
tary cross section. As the elementary cross section, we 
have used the experimental data of p + d ^ t + tt^ re- 
action reported in Ref [18] in the CM frame by isospin 
symmetry. We parameterize the data shown in Fig. 5 in 
Ref [H] as. 



CM 



da\ 
^Aic 



^0-0.852 cos(.-e--i)+0.051 [^^/sr] , (14) 



and transfer the cross section to that in the lab frame by 
Eq. ([S]). We show the values of elementary cross sections 
used in this paper in Table [TJ 

We, then, show the differences of the kinematics be- 
tween the elementary process d + n ^He + vr^ and the 
meson bound state production A{d, ■^He)(yl— 1)^tt~ with 
a heavy nucleus target. We show in Fig. [1] the momentum 
transfer q at O^^He ~ ^° ^^'^ 2° for the elementary process 
and the reaction with the heavy nuclear target for pion 
production with the meson binding energy B.E. — and 
the neutron separation energy Sn = 0- We find that the 
recoilless condition (|q| = 0) is only satisfied at 6]^^^ = 0° 
at Td — 455 MeV for both reactions. At different incident 
energy, the momentum transfer takes obviously different 
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Fig. 1. Momentum transfer q (upper) and kinematical cor- 
rection factors K (lower) of the (d^Yie) reaction at ^Jj^^, — 0° 
(thick lines) and 2° (thin lines) plotted as functions of the inci- 
dent deuteron kinetic energy T^. In the upper frame, the solid 
lines show the results for the formation of tt meson in the heavy 
target nucleus, and the dashed lines those for the elementary 
process. The pion binding energy B.E. and the neutron sepa- 
ration energy S„ from the target nucleus are assumed to be 
MeV. The vertical line indicates the incident energy satis- 
fying the recoilless condition at ^Jj^^, = 0°. 
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Fig. 2. Momentum transfer q (upper) and kinematical correc- 
tion factors K (lower) of the (d,^He) reaction at Q"^^^^ = 0° and 
2° for the formation of tt meson bound states in the heavy tar- 
get nucleus plotted as functions of the reaction Q-value. The 
incident deuteron kinetic energy is fixed to be Td = 500 MeV. 
The neutron separation energy S„ is fixed to be 8.8 MeV (solid 
lines) and 10.2 MeV (dash-dotted lines). The vertical lines in- 
dicate the TT meson production threshold for Sn = 8.8 and 10.2 
MeV cases. 



value for these two reactions. We also show the incident 
deuteron energy dependence of the K factor in Fig. [TJ 
We found that the K factor is 1 for the recoilless kine- 
matics with B.E. = and Sn = as naturally expected. 
The K factor is at = 407 MeV and = 430 MeV 
for 6'Jfi5?^ = 0° and 2°, respectively. These energies cor- 
respond to the pion production threshold for the elemen- 
tary process. And the K factor increases with the incident 
deuteron kinetic energy T^. 

In Fig. [21 we show the momentum transfer for the nu- 
clear target case and the K factor for = 500 MeV at 
^dHe ~ ^"^"^ ^° functions of the reaction Q-value, 
which is the kinematics considered in the missing mass 
spectroscopy of the deeply bound pionic atoms in this ar- 
ticle. We find that the K factor gradually varies about 
10% within the Q-value range considered here for both 
^dHc ~ ^° ^^"^ 2° cases, and the K factor decreases about 
20% at e^l^^ = 2° from the value at 6^1^^ = 0°. Thus, 
we think the K factor should be introduced to study the 
angular dependence of the (d,'^He) spectra. 



In Fig. [21 we show the calculated (d,'^ He) spectra at 
finite angles for the bound pionic states formation. Two 
kinds of the neutron wavefunction have been used for the 
calculations, which are the harmonic oscillator wavefunc- 
tion and the calculated wavefunction using the neutron 
potential reported in Ref . [IQi . Though the absolute value 
of the calculated cross sections depend on the neutron 
wavefunction used as shown in Ref. |12) . the angular de- 
pendence of the spectra resembles each other. We find 
that the spectra have a strong angular dependence and 
the shape of the spectra are much different at finite angles 
from that at 0°. The largest peak structure at Q = —137.8 
MeV in the forward spectra is strongly suppressed at finite 
angles and the spectra show the structure of three peaks 
at 0^^^^ > 2°. The overall strength of the spectra has also 
angular dependence and is smaller for larger angles. 

In Fig. [31 we show the dominant subcomponents of 
(dj^He) spectra for each scattering angle O]^^^. Here, we 
have used the neutron wavefunction obtained by the po- 
tential in Ref. [H] and we have also shown the spectra 
with better energy resolution 100 keV case. At O]^^^ ~ 
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Fig. 3. Calculated ^^^Sn(d,'^He) spectra for the formation of 
the pionic bound states at 6^^^^ = 0° (solid lines), l°(dash 
lines), 2° (dash-dotted lines) and 3° (dotted lines) plotted as 
functions of the reaction Q-value. The harmonic oscillator (up- 
per) and the phenomenological [15] (lower) neutron wavefunc- 
tions are used. The incident deuteron kinetic energy is fixed 
to be Tfj = 500 MeV. The instrumental energy resolution is 
assumed to be 300 keV FWHM. The vertical line indicates the 
pion production threshold Q = —141.6 MeV. 



0°, since the reaction is close to recoilless, the peaks of 
[(Is)^(g) (3si/2)~^] and [(2s)^ (g) (3si/2)^^] subcomponents 
appear clearly in the spectra for both energy resolution 
cases. The contributions of (3s)^ and (4s)7r also show the 
isolated peaks for better energy resolution (100 keV) case 
even though they are small. At (i^^^^ = 1°, the contribu- 
tion for [(2s)7r ® (3si/2)^^] is suppressed and can only be 
seen clearly in a better energy resolution case. At larger 
angles, the pionic (2p) state contributions become rela- 
tively larger and dominate the peak structure around Q = 
-139.4 MeV and -140.4 MeV. We can expect to observe 
the peak structure composed from [(2p)7r ® (3si/2)^^], 
[(2p)jr ® (2^3/2)"^] and [(2p)^ ® (^h^\li)n^\ subcompo- 
nents. Though, the separation energies of these 3 neutron 
levels differ from each other only within 60 keV [T2] and 
their contributions can not be distinguished, we can ex- 
pect to deduce the information on the pionic 2p state. 
These contributions of pionic 2p state can not be seen in 
the spectra at B"^^^^ = 0° since they are hidden in the tail 
of the large Is contributions. At finite angles, due to the 
significant suppression of the Is contributions, the 2p con- 
tributions can be observed even if they are a little smaller 



than those at the forward angle. Thus, to observe the spec- 
tra at finite angle is valuable. 

Finally, we comment on the small peak structures in 
the pionic unbound energy region which appear in bet- 
ter energy resolution cases in Fig. |4] right panels. They 
are the contributions of lightly bound pionic atom forma- 
tion with one of the deep neutron hole states correspond- 
ing to the excited states of the daughter nucleus such as 
[(3p)^0(2d5/2)-i], [(4p).(»(2d5/2),T'], [(3d).<»(2d5/2)-i], 
[(4c?)7r ® (2^5/2)^^] and so on. The wavefunctions of the 
lightly bound pionic atom have large spatial dimensions 
and small overlaps with nucleon wavefunction, and thus, 
both widths and formation cross sections of these bound 
states are small. Therefore, we can observe contributions 
of these states formations only in the spectra with better 
energy resolution. 



4 Summary 

We study the formation of deeply bound pionic atoms 
in the (d,'^He) reactions theoretically and show the an- 
gular dependence of the expected spectra, which can be 
observed in experiments. 

We develop the formula to include the kinematical cor- 
rection factors to the effective number approach to obtain 
more realistic angular dependence of the (d,'^He) spectra. 
We show the behavior of the kinematical correction factor 
and the angular dependence of the -'^^^Sn(d,'^He) spectra 
at Td = 500 MeV for the formation of the pionic atoms. 

We find that the spectra are dominated by the subcom- 
ponents including (2p)^ state at larger scattering angles 
— 2°, while they are dominated by the (Is)^ and 
(2s)7r states at forward angles. The peaks are well isolated 
and can be observed in the experiments with the good en- 
ergy resolution. Thus, we can conclude that we can obtain 
information of deeply bound pionic 2p state in addition to 
Is and 2s states by observing the spectra at finite an- 
gles. As indicated in Ref. [12], the observation of several 
deeply pionic bound states in a certain nucleus will help 
to deduce precise information of pion properties and the 
chiral dynamics at finite density [6l[7l|8]. We believe that 
our results provide a good evaluation for further experi- 
mental studies of the states reported here, which should 
contribute to the development of the field. 
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Fig. 4. Calculated ^^^Sn(ci,^He) spectra for the formation of the pionic bound states at 6]iiic = 0°, 1°, 2° and 3° are plotted as 
functions of the reaction Q-value. Dominant subcomponents [{nl)^ (gi (n£j)~^] are indicated in the figure. The phenomenological 
neutron wavefunctions [19] are used. The instrumental energy resolution is assumed to be 300 keV FWHM (left) and 100 keV 
FWHM (right). The vertical line indicates the pion production threshold Q — —141.6 MeV. 



